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Abstract. In this note we provide dispersive estimates for Fourier integrals 
with parameter-dependent phase functions in terms of geometric quantities 
of associated families of Fresnel surfaces. The results are based on a multi- 
dimensional van der Corput lemma due to the first author. 

Applications to dispersive estimates for hyperbolic systems and scalar higher 
order hyperbolic equations are also discussed. 

1. Introduction. Dispersive estimates for solutions to linear evolution equations 
are vitally important for the study of a wide range of problems, including ques- 
tions of stability of solutions to corresponding non-linear equations, descriptions of 
large-time asymptotic profiles of solutions or scattering results based on the time- 
dependent approach. 

For the particular situation of hyperbolic evolution equations with a hyperbolic 
time-asymptotic profile, such estimates can be based on corresponding decay es- 
timates for Fourier integral operators with real phases and certain well-controled 
amplitude functions. The aim of this note is to provide some robust estimates for 
such Fourier integrals in terms of geometric properties of associated Fresnel surfaces. 

2. Fourier integral operators, notation and main results. We are interested 
in L^-L 00 estimates for a family of Fourier integral operators depending on a pa- 
rameter t £ R+ (thought of as time later on) 

T t u(x)= [ e^< +t ^ x ^a(t,x,0u(m, (1) 

JR™ 

where a(t, £,£) is uniformly in t a (global) symbol of order zero, 

|^o(t,x,OI<C a <0~ H . W 
with (£) = yl + |£| 2 , and <p(t, x, £) is a homogeneous real- valued function satisfying 

C- 1 |£|<0(t,x,O<C|£|, \d^(t,x,0\ <C' a \e- H (3) 

uniform in t and x. For simplicity, we always assume u £ 5(R n ), estimates extend 
by continuity arguments to larger spaces. If we are only interested in estimates of 
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Ttu in the L°°-norm, we can treat t and x both as parameters and derive estimates 
for the corresponding stationary model operator 

T:u^ f e i( ^ +<w?)) a(£HOde (4) 



A precise understanding of the appearing constants in the estimate translates to 
estimates of the t-behaviour of the full model (1). 

2.1. Fresnel surfaces and contact indices. We turn to the analysis of (4) and 
define the geometric quantities which determine its L°° behaviour. The consider- 
ation follows essentially Sugimoto [15], [16], for the non-degenerate situations see 
the earlier works of Strichartz [14] or Brenner [1] . 

We assume that the phase <^>(£) is non- negative, 1-homogeneous and smooth, and 
define the associated Fresnel surface (also called slowness surface) 

E = {£ £ R" : 0(C) = 1}. (5) 

We introduce two indices for such a Fresnel surface E, assuming that it is of class C k 
with k being sufficiently large. For p £ E we denote by T p the tangent hyperplane 
to E at p. Then for any plane H of dimension 2 which contains p and the normal of 
E at p we denote by 7(£,p, H) the order of contact between the curve E n H and 
its tangent H n T p . Furthermore, we set 

7(E) = sup sup j(E,p,H), 70(E) = supinf 7(E,p,i7). (6) 

p H p H 

The definition implies directly that 2 < 70(E) < 7(E). Furthermore, for isotropic 
problems E is a dilation of the sphere S n_1 and 7(§™~ 1 ) = 7o(§"~ 1 ) = 2. Moreover, 
if the Gaussian curvature of E never vanishes E is convex and we have 7(E) = 2. 

These indices arc directly related to mapping properties of (4). The following 
theorem is a combination of results from [15], [16] in the analytic setting, improved 
by the author [9], [10] to the smooth setting as well as to the limited regularity 
cases. Moreover, the corresponding results in [9], [10] do not require the positive 
homogeneity of the phase in £ and allow the dependence on parameters. 

As usual, Bp 9 (R") denotes the Besov space of regularity r modelled over LP . 

Theorem 2.1. Assume that <\> is smooth on M. 1+n x (W 1 \ {0}), positive and 1- 
homogeneous in £, E smooth and a £ S® (M. n ). Then the operator T defined by (4) 
maps 

B[ )2 (K™) L°°(R n ) (7) 
where r = n — if E is convex and r = n — 7o ( S ) if E is arbitrary. 

The operator norm in this estimate depends on finitely many symbol seminorms 
of a, the estimates of 4> and its derivatives together with some quantitative measures 
of the contact between E and tangent lines. 

To control the order of contact of E by tangent lines quantitatively, we will 
introduce another quantity x(E). First, assume that E is convex and that it is 
of class (7 7 ' E ) +1 . For p <E E, rotating E if necessary, we may assume that it is 
parameterised by points {(y,h(y)),y G 0} near p for an open set 17 C R n_1 . For 
p = (y,h(y)), let us define 



dpi 



h(y + pu)\ p=0 



(8) 
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From the definition of 7(E) it follows that xr(E,p) > for all p 6 f2. Indeed, from 
the definition of 7(E,p, -ff) it follows that if ui is such that y + put e H, then 

a7(s,p,-ff) 

x(X,p,H) = g T(e<Rfl) % + H[H >0- (9) 

Now, we clearly have Tnph(y + pw)| p =o > *c(E,p, H), and hence we have 

x(E,p) > since the set |w| = 1 is compact. Noticing that >r(E,p) is a continuous 
function of p, by compactness of E it follows that if we define 

x(E) = minxr(E,p), (10) 
pes 

then >f(E) > 0. If E is not convex, we define similarly 



j<o(E) = min sup 



pes 



70 (S) 

E 

J=2 



g-h(y + pco)\ p=0 



(11) 



Again, we have >ro(E) > 0. 

Considering again the example of spheres E = rS n_1 of radius r, we see that 
7(E) = 70(E) = 2 and that >r(E) = >fo(E) ~ - is a measure of scalar curvature. 

2.2. Main results. Estimates of (1) depend on indices of a family of Fresnel sur- 
faces defined in terms </>(t, x, £), 

E t , x = U e R" : 0(t, = 1}, t > to, x£l". (12) 

For a first result we assume that E^ is convex for all t > to and all x €E K™. 

Theorem 2.2. Assume that cj>(t,x,£) is real-valued, continuous in t and x and 
smooth in £ S R n \ {0}, homogeneous of order one in £ and suc/i t/iat /or some 
to > and C > we have 

C- 1 ^ <</>(!, x,t)<C\t\ and |^(t,x,OI<C|f| 1 -H (13) 

/or a/Z t > to, aZZ £ 7^ and all multi-indices a with \a\ < max{7 + 1, |_(rz — 
1)/tJ + 2}. Assume further that all the sets T^t,x defined by (12) are convex for all 
t > to and assume that 

limsupsup7(Ei >x ) < 7 together with liminf inf x(Ej )X ) > 0. (14) 

Suppose further that the amplitude a(i, x, £) satisfies 

\d£a(t,x,0\ < Ca(£)~ W /oraZZ |a| < L(n - 1)/tJ + 1. (15) 
T/ien the operator family T t defined by (I) satisfies for all t > to the estimate 

\\Ttu\\ Lx ^) <Ct ^HuHBr^n) (16) 

with r = n — — - . 

7 

In the non-convex situation decay rates can be much lower. Results can be 
improved by using a more detailed analysis of the geometric situation near points 
of high contact order. 

Theorem 2.3. Assume that </>(t,x,£) is real-valued, continuous in t and x and 
smooth in £ € K n \ {0}, homogeneous of order one in £ and such that for some 
to > and C > we have 

C-^\ <4>(t,x,0 <C|£| and \d^(t,x,0\<C\e^ al (17) 
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for all t > to, all x, £ ^ and all multi-indices a with \ot\ < 70 + 1- Assume further 
that all the sets Ti t , x defined by (12) satisfy 

limsupsup7o(E tia ;) < 70 together with liminf inf >«b(Ej iX ) > 0. (18) 

t— >oo x i— >oo X 

Suppose further that the amplitude a(t,x,£) satisfies 

\8£a(t,x,0\<C a (Z)-M for all \a\ < 1. (19) 
Then the operator family T t defined by (I) satisfies for all t > to the estimate 

llrtltlU-^n) < Ct~™ ||u|| fl r ^n) (20) 

with r = n — — . 

7o 

The proofs of both theorems are similar to the ones in [13] and full details are 
omitted here. The dependence on x is parametric, and the assumption assure that 
the resulting estimates are uniform in x. 

If we make further regularity assumptions with respect to the a;- variable it follows 
that the operators T t are uniformly i 2 -bounded such that interpolation implies 
the usual L p -L q decay estimates. For completeness, we include a corresponding 
statement. The assumptions are chosen such that after substituting x = ty and 
t£ = r\ Theorem 2.5 of [11] applies: 

Theorem 2.4 ([11]). LetY y ,Y^ C W 1 be open cones. Let operator T be defined by 

Tu(x)= [ [ e*(** + *<»>e»a(y,t)u(3,) dy d£, (21) 

where a(y,£) £ C°°{W> x K£), supp a C T y x T i; and <p(y,Q E C°°{T V x T e ) is a 
real-valued function. Assume that 

\d^a(y,0\<C aP , 

for \a\, |/3| < 2n + 1. Also assume that 

\detd y d£<p(y,£)\ > C > on T y x r 6 

that 

\dt<p(,x,Z)-d e <p(v,S)\>C\x-y\ for x, y € T y , £ e T e , 

and that 

\dyd£<p(y,€) <C a \, \d y dl<p(y,£)\ < C p on supp a 
for 1 < |a|, \/3\ < 2n + 2. Then the operator T is L 2 (W n ) -bounded, and satisfies 

||r|| L 2^ L 2 < C sup Wi |0|< 2n+1 ||d^a(y,Ollz^(R£xMp- 

We can now formulate the corresponding L 2 -result: 

Theorem 2.5. Let 

T t u{x)= f ^*-W>**»a(t,z,t)G{tte, ( 22 ) 
Assume <f>(t,x,£) is real-valued, l-homogeneous in £ and satisfies 

\det(I + td x ds<f>(t,x,0)\ > Co >0, \t^d^d^(t,x,0\ < C aJi (23) 

uniform in t > to, x, £ ^ and for 1 < \a\, \/3\ < 2n + 2. Assume further that 
a(t,x,t;) is supported in t|£| > C for some constant C and that 

\tWdZdj!a{t,x,0\<C a ,p (24) 
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for all \a\, \/3\ <2n + 2 and t > to. Then the family Tt in (22), t > to, is uniformly 
L 2 -bounded, 

\\T t u\\ L 2 {Rn) <C\\u\\ L 2 {Rn) . (25) 

Let us sketch how Theorem 2.4 implies Theorem 2.5. Using homogeneity, we can 
rewrite (22) as 

T t u(x)=t-" [ e i ( z -' )+<# ' t (^^ , '»a t (t,z,?7)2(r//i)d?7, (26) 

JR™ 

with x = tz,t^ = 77, and 

4>t (t, z, rf) = 4>(t, tz, n), a t (t, z, n) = a(t, tz, n/t). 

By rescaling, the powers of t in the L 2 -estimate cancel, and it can be shown that 
the adjoint to the operator in (26) satisfies assumptions on Theorem 2.4 uniformly 
for t > t . 

We also note that the assumptions on the mixed derivatives of 4> as m (23) can 
be relaxed by looking only at |a| = 1 or (3 = 1, as in Theorem 2.4. 

2.3. Appendix: Fourier transform of surface carried measures. Let £ C 

M. n be a closed hypersurface. Then the indices 7(S) and 70 (£) arc related to 
decay estimates for Fourier transforms of surface carried measures. Consider for 
/ G C°°(S) the integral 

u(x) = f e- «/(Ode, (27) 

where d£ denotes the (n— l)-dimcnsional surface measure on E. Then the following 
theorem is valid. 

Theorem 2.6. 1. Assume E is a smooth hypersurface and convex. Then 

\u(x)\<C(xy^m\\f\\ Crm (28) 

with r>(n- 1)/7(E) + 1. 
2. For general (non-convex) smooth hypersurface £ it follows that 

\u{x)\<C{x)-^\\f\\ c ^y (29) 

The proof of the first statement in [15] uses real analyticity of the surface, the 
stronger results are due to methods of [9] and [10]. Smoothness of the surface can 
be further relaxed to the class C 7 with 7 = 7(E) (or 70(E), respectively). 

3. Applications to hyperbolic evolution equations. The estimates discussed 
so far have applications for dispersive estimates of hyperbolic evolution equations. 
Without aiming for completeness we will discuss some possible applications. 

3.1. Wave models. As a first example we consider wave models with time-dependent 
coefficients, 

n 

d\u- ^ a tj (t)d Xt d Xj u = (30) 

in R x R" and with Cauchy data u(0, ■) = uq, u t (0, ■) = u\, both assumed to be 
Schwartz for simplicity. Denoting 

T{£} = {ae C°°(R) : \d?a(t)\ < c fc (t) £ - fc }, (31) 
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we assume a%j G 7~{0}, a,ij(t) = a,ji(t), together with the strict hyperbolicity as- 
sumption 

n 

a(t,$='£a ij (t)S i Z j >C\t\ 3 , (32) 

ij = l 

for some C > uniform in t. Within the zone {(1 + > 1} solutions to (30) are 
given by Fourier integrals of the form (1) with phase functions 

4>±{t,£) = ±\ I yffiffijM. (33) 
1 Jo 

In this situation the Fresnel surfaces are given by 

E t = U G K" : t = J v^Od^j (34) 

and satisfy 7(St) = 2 together with uniform upper and lower bounds for x(£t). 
This follows directly from [6, Lemma 3.1]. Thus Theorem 2.2 is applicable and in 
combination with a uniform bound on the £ 2 -energy one concludes that solutions 
to (30) satisfy 

\\d t u(t, ■% + ||Vu(t, -)||, < C pq (l + i)-^(?-i) (\\ Ul \\ WP , rp+ i + ||u 2 |kp.rp) (35) 

with l<p<2<q< oo, pq = p + q and r p = n(l/p — l/q). 

For a detailed consideration of such models we refer to Rcissig [6], [7]. It is 
possible to include further lower order terms b(t)dtU, J2j c j{t)d Xj u and d(t)u with 
b,c.j G T{ — 1} and d G T{— 2} (together with further assumptions to control the 
low-frequency behaviour). We refer to [17] for an example of such a treatment. If 
lower order terms belong to different classes, the asymptotic behaviour of solutions 
might change its type, see, e.g., [18] a discussion about dissipative wave equations 
and effective dissipation. 

The assumptions on coefficients can be weakened to classes 

%{£} = {a G C°°(R) : \d?a(t)\ < c k (t log(e + t) v ) l ~ k } (36) 

with parameter v G [0, 1] where v G (0, 1) implies a small loss of decay and v = 1 
yields a polynomial loss. It can be shown that without further assumptions v > 1 
leads generically to supra-polynomial growth of the energy. Additional assumptions 
to compensate this behaviour for wave equations with variable speed of propagation 
were considered by Hirosawa [3], see also [19], but the estimates of Theorems 2.2 
and 2.3 need to be refined in order to obtain dispersive estimates in this situation. 

3.2. Higher order scalar equations. In a similar way we can treat homogeneous 
higher order equations of the form 

m— 1 

D ™ U +5Z J2 afc,«(*) D «D> = 0, (37) 

k=0 \a\=m—k 

Diu(0,-)=Uj, i = 0,l,...,m-l, (38) 
where D = —id, a,k.a G T{0} and the characteristic equation 

m— 1 

Tm +E E ^At)r k C = (39) 

k—0 \a\— m — k 
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has m uniformly distinct real solutions n(t, £),. .., T m (t, £), \n(t, £) — Tj (t, £) > C|£ 
for i 7^ j and some C > 0. Under the stronger assumption that a' k a (t) G L 1 (R+), 
this model was treated by Matsuyama and the first author, [5], the full form follows 
from [13]. If we assume 



E 



(16* 



< C (40) 



uniform in j, t and £ ^ 0, solutions are given as sums of Fourier integrals of the 
form (1) with phases 

1 



<f>j(t,t)=tJ-rjV,t)M (41) 

and uniformly bounded amplitudes. Again, the results from Section 2.2 provide 
dispersive estimates for solutions to such problems. We remark that for highcr 
order equations the associated Fresnel surfaces T, t need not be convex and therefore 
the full generality of Theorems 2.2 and 2.3 is of importance here. 

We note that the analysis in [5] is based on the asymptotic integration method 
for constructing solutions for equations with homogeneous symbols, and it requires 
only C^-regularity of coefficients, i.e. only that ak, a £ C 1 in terms of the regularity 
in time, which is crucial in applications to the Kirchhoff equations, see e.g. [4]. 

3.3. Hyperbolic systems. In [12], the authors announced results on dispersive 
estimates for a general class of pscudo-diffcrcntiablc hyperbolic systems with time- 
dependent coefficients and arbitrary 7. D'Abbico, Lucente and Taglialatcla, [2], 
considered differential hyperbolic systems with 7 = 2, while in [13] the authors 
gave a rigorous treatment for general pseudo-diffcrcntiable hyperbolic systems with 
time-dependent coefficients and arbitrary 7. For the sake of simplicity we formulate 
the results for the homogeneous situation. We consider 

T> t U = A(t,D x )U, U(0 r ) = U (42) 

for a m x m matrix A(t,T) x ) of f-dependent Fourier multipliers. We denote 

5{m 1 ,m 2 } = {o(t,0 £ C°°(IR x (R n \ {0})) : \^- m ^a(t, £ TW}, (43) 

the condition being uniform in £. We assume that the symbol Aft,^) € 6>{1,0} is 
1-homogencous in £, a self-adjoint matrix and its eigenvalues Ti (£,£),. . . , T m {t, £) 
are uniformly distinct, \n(t, £) — Tj(t, £)| > C|£| for i ^ j and some C > 0. Without 
further assumptions it follows that the solutions to (42) can be represented in the 
form (1) with phases given by (41). 

Assumptions can be relaxed, without assuming self-adjointness we have to require 
real eigenvalues together with a uniform symmetrisability of the matrix A(t, £) to- 
gether with further assumptions to control the large-time influence of lower order 
terms. We refer to [13] for details. 
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